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Abstract. We study by theoretical analysis and by direct numerical simulation the dynamics of a wide 
class of asynchronous stochastic systems composed of many autocatalytic degrees of freedom. We describe 
the generic emergence of truncated power laws in the size distribution of their individual elements. The 
exponents a of these power laws are time independent and depend only on the way the elements with 
very small values are treated. These truncated power laws determine the collective time evolution of the 
system. In particular the global stochastic fluctuations of the system differ from the normal Gaussian 
noise according to the time and size scales at which these fluctuations are considered. We describe the 
ranges in which these fluctuations are parameterized respectively by: the Levy regime a < 2, the power 
law decay with large exponent (a > 2), and the exponential decay. Finally we relate these results to the 
large exponent power laws found in the actual behavior of the stock markets and to the exponential cut-off 
detected in certain recent measurement. 

PACS. 05.40.+j Fluctuation phenomena, random processes, noise, and Brownian motion - 05.70.Ln 
Nonequilibrium thermodynamics, irreversible processes - 02.50.-r Probability theory, stochastic processes, 
and statistics 



1 Introduction 

It was realized since a very long time that the fluctua- 
tions r(f) of stochastic systems made of many degrees of 
freedom are not generically distributed by Gaussian prob- 
ability distributions G|,0|. On the theoretical side, P. Levy 
H discovered at the beginning of this century that the 
central limit theorem allows for a family of distributions 
which decay at infinity as a power law: 



P(r) 



with < a < 2 



(1) 



These distributions can be thought as the limit distribu- 
tions for random walks with steps of sizes Sj distributed 
by a power law 

P(s) ~ s- 1 -". 

Such processes were named Levy flights pj . 

In nature, the situation turns out to be more compli- 
cated: many of the measurements produced probability 
distribution functions which look like Levy distributions 
for a certain range of the stochastic variable but are cut 
off. That is, they change their character above a certain 
threshold, e.g., by becoming exponentials or changing to a 
power decay with a > 2 as found for the returns of stock 
market data Particular examples of such distri- 

butions were studied under the name of truncated Levy 
flight ||. As required by the central limit theorem, for 



very large time intervals t, these distributions cross-over 
into a quite Gaussian behavior. 

Recently, both the power law probability distribution 
of the individual steps and the truncated Levy distribu- 
tion of the fluctuations were explained generically by a 
series of Generalized Lotka-Volterra (GLV) models P,|lC|, 
|ll],[l2). These models represent realistic financial, biolog- 
ical and social systems composed of many autocatalytic 
and competing stochastic subsystems [lC( ]. 

In this paper we study the exact nature of the cut- 
off region and the details of the cross-over process in the 
framework of the GLV models. We especially describe the 
emergence of the tail distribution with a « 3 power law 
and the exponential cut-off. The theoretical analysis below 
yields well defined quantitative predictions on the various 
temporal and geometric properties of the probability dis- 
tribution functions, which are verified by the numerical 
simulations and compared with the actual measurements 
of the Hong Kong stock market. Most of the results extend 
however to other models too. 



2 Theoretical analysis 

2.1 Autocatalytic systems and power law 

The analysis here applies to a large range of dynamical 
models JL0| . For definiteness we consider a specific system 
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W{t) formed of subsystems w^t), i = 1,2, 2V: 

W(t) = tfli(t)+ioa(t) + ... + «;*(*) (2) 

distributed by a power law cumulative distribution func- 
tion: 

Pi>w) = [w/w min (t)}- a (3) 

with the corresponding density P(w) oc u> _1 ~ Q , where 
Wmin is some lower cut-off 



Wmin(t) = gw(t), 



(4) 



i.e., usually given in terms of a fraction q of the average 
tD(i) = W{t)/N. 

Such distributions emerge naturally in autocatalytic 
systems |)|,[l0],[ll],|l2j of the generic form: 

Wi(t + 1) = X(t)wiit) + aw{t) - bwi(t)w(t). (5) 

For the purpose of this paper we describe them by a simple 
dynamics P,p2[ which consists in choosing randomly at 
each time stepfor updating one of the Wi 's and multiplying 
it by a random factor A(t) extracted at each time from a 
probability distribution p(A): 



Wiit+l) = \(t)Wiit), 

with the lower cut-off 

Wiit + 1) > qw(t). 



(6) 



(7) 



That is, the updated variable Wi(t + 1) is constrained to 
be not smaller than a lower bound qw(t), i.e., 

Wiit + 1) = max{iUi(i + 1), gtD(i)}, 

where w(t) is the average value calculated at earlier time 
t. Note that all the w^s are strictly positive and therefore 



Wiit) < W(t) = Nw(t). 



(8) 



In systems of the type Eqs. (||) and (0), it has been 
shown y,|l2) that even if the system is not stationary the 
power law Eq. ((H) holds, and for given N and q in the 
range 1 > q 1/lniV, the exponent a is given by the 
relation: 

a -1/(1 -q) (9) 

Usually in financial applications q ~ 1/3, and then the 
exponent a ~ 1.5 within the stable Levy regime. 

During our analysis of this section we will take the 
distribution /c(A) centered around 1: 

X(t) = l+git), (10) 

and consisting of just 2 equally probable values 

git) = ±M. (11) 

The conclusions of our analysis are not changed if one 
assumes an arbitrary distribution p(A) with finite strictly 
positive support. 



Here we are interested in the distribution of the "re- 
turns" : 

r(r) = [W(t + T)-W(t)]/W(t) (12) 

as a function of the time interval r. The use of the term 
"return" in ( |T^ ) as a measure of the fluctuation in the 
system is borrowed from the financial applications where 
Wit) is the stock market index and r(r) is the relative 
gain/loss one incurs after a time lapse r. 

Since the variation of Wit) at each time coincides with 
the variation of the individual Wi(t) which happened to 
be randomly selected for updating by Eq. (||) at time t, 
the value of r(r) as defined in Eq. ( |l2|) is the result of a 
random walk 

(13) 



r(r) =5~J«i(t + fr), 



with steps Siit + k), k = 0, 1, r — 1 of sizes ]l3| : 

Si(t + k) = [wiit + k + l)- w^t + k)]/W{t), (14) 

which according to Eqs. (||), ([iTj| ) and ( [j~L| ) is written as: 

Si(t + k) = ±Mwi(t + k)/W(t). (15) 

Therefore, the sizes of the (absolute values of the) individ- 
ual steps Si(t + k) in the random walk/flight process r(r) 
(Eq. (|13|) ) have a probability distribution function similar 
to Eq!"7|) (up to the factor M/W): 

Si> s) = [s/s min )- a , (16) 

where s min = Mw min /W = Mq/N. 

2.2 Origin of truncation in autocatalytic systems 

A crucial fact of the present paper is that the power law 
distribution Eq. ( |l6| ) from which the individual steps Eq . 
( |l4|) composing the random walk r(r) through Eq. (fLq ) 
are selected is truncated from above |14[| . 

Indeed, since Wi < W for any i as shown in Eq. (^|), 
the individual steps Eq. (|l4|) cannot be larger than the 
fixed value M, i.e. 



S(> M) = 0. 



(17) 



Note that this bound in the size of the individual steps 
of the random walk r(r) does not depend on the number 
N of participants in the game nor on the exponent a of 
the power law, nor on the lower cut-off w m i n . One cannot 
therefore hope that the effects of this upper bound would 
somewhat become irrelevant, and as seen below (Eq. (p6|)), 
by increasing N one can only delay the time 

1 max 1 y 

by which these effects become dominant. Therefore, for fi- 
nite time intervals r and infinite N, the effect of the trun- 
cation disappears (together with the vanishing of the am- 
plitude of elementary fluctuations |r(l)| < N 1_a JTo| ) . In 
finance, however, we are generically in the opposite limit, 
e.g. a finite number of traders trading for very long time 
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intervals r which allow them to perform a total number 
of elementary trading operations much larger than T max . 

Our goal in the sequel will be therefore to follow in 
detail the process by which the r = 1 (truncated at r > 
M) power distribution 

R(> r, 1) = S(> r) = [r/s min ]- a = [r/(qM/N)]- a (18) 

evolves for increasing time interval r, and then approaches 
the infinite time (r 3> T max ) Gaussian distribution. We 
will analyze explicitly R(> r, r) only for positive r's to 
avoid unnecessary complication. However, the analysis for 
negative r values is very similar. 

In fact, for the probability distribution function of re- 
turns R(> r, t) one obtains a symmetric probability den- 
sity R(r,r) = dR(> r,r)/dr which for r = 1 coincides 
with the S(s) = dS(> s)/ds probability density. On a 
log-log scale the probability density lni?(r, 1) vs. lnr is 
a line ending sharply around r = s m in- We will see that 
for larger r values the sharp tip will erode into a flatter 
"dome" and the complete/exact vanishing of the distribu- 
tion at the upper cut-off r — M will evolve into a steep 
but continuous decay. 



2.3 Power law and truncation for r < r max 

The probability R(> r, t) to reach after a r-steps walk a 
distance r or larger is of course a result of the probabil- 
ities of the r individual steps which compose the "walk- 
ing" /" traveling" /" flying" p process. Therefore its char- 
acteristics depend on the 3 crucial properties of the indi- 
vidual steps distribution S(> s) obtained from Eqs. (jig)- 

&■ 

1. The great majority of the individual steps Si are of 
order s mean (the average of Si) and less, i.e., in the 
range: 

Smin = qSmean = qM / N < S < 2s mean = 2M/N. 

2. The steps of larger sizes (say larger than 2s mean ) are 
very rare: 

R{> 2M/N, 1) - (q/2) a < 0.07 (20) 

(We take a ~ 1.5 as in real wealth distribution [p^[). 

3. There are no individual steps in r (r) of size larger than 
M: 

R(>M, 1) = 0. (21) 

Due to those properties, as one increases the time inter- 
val r from 1, the initial (truncated) power law distribution 
Eq. ( p^[ ) R(> r, 1) = S(> r) is not significantly affected 
for most of the r range. For small r's, the corrections to 
the power law are in fact limited only to the lower and 
upper cut-off regions and are analyzed below. 



2.3.1 The low r "dome" -like region 

The low r region in i?(> r, r) is affected even for small 
r > 1 because there is a large probability that all of the r 
steps are of the order s mean and lower. Consequently, the 
probability R(> r, t) for values r < Ts rnean — tM/N is 
not given anymore by the probability R{> r, 1) = 5*(> r) 
of obtaining it through a single step, but rather by a com- 
binatoric sum of probabilities of having r (small, pos- 
itive and negative) steps summing up to r. This is of 
course very similar to the way one estimates (through 
Poisson/Binomial expansion) the probability of a distance 
r after a r steps Gaussian walk. The consequence is a 
smoothening of the sharp end at R(r = 0, r). This con- 
cretizes in the appearance of a "dome" -like shape in the 
central region (around r = 0) of the R(r, r) probability 
density. 

To estimate (the time dependence of) the extent of 
the "dome" -like region, we demand that the probability 
of achieving distances r through 2 or more steps is larger 
than the probability of achieving it through one step. Since 
(rS(r/2)) 2 is the probability of having during r steps at 
least 2 steps of sizes at least r/2, and rS(r) is the proba- 
bility of having at least one of the r steps of size at least 
r, the condition describing the "dome" -like region is: 

(rS(r/2)) 2 > rS(r), 

which by substituting the power law Eq. (|l6| ) for S(> s) 
becomes: 

{T[r/{2s mm )]- a } 2 > T(r/s mrn )- a , 

gives the central dome region: 

r < 4s min r 1/Q = AqM/Nr 1 ^. (22) 

since q is in (finance markets) practice not far from 1/4 
we denote 

Tdome = S mean T X l a = M/Nt 1 '" (23) 

and with this notation, the condition Eq. ( p2| ) becomes 

r < Tdome- (24) 

Thus, the power law Eq. ( |l8| ) remains unchanged in the 
range 

Tdome <r < M, (25) 

and the power law region disappears completely when the 
upper and lower limits of this interval coincide: 

Tdome 

i.e. (cf Eq. ( ^3|) by the time that 

r 1/a M/N = M 

holds. This gives the maximal time for which one still has 
a power-like region: 

Tmax = N a . (26) 
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Since the power law region is the crucial feature of the 
(truncated) Levy distribution, T max is essentially the max- 
imal time interval for which the returns ( |l2[ ) still maintain 
a Levy-like probability distribution. Note that this value 
for T max agrees with the early estimations of Ref. B based 
on the scaling i?(0, t) ~ t 1 /™ of the probability distribu- 
tion peak. 



2.3.2 The extremely large r > M region 

In the upper cut-off region r > M the probability distri- 
bution function R(> r, r) is affected by increasing t from 
1 to larger values. Indeed, while values r > M are com- 
pletely disallowed for r = 1 because of the truncation Eq. 
(pl|), for t > 1 one can have r's in the range M < r < tM. 
The probability of obtaining such values of r corresponds 
to the probability of selecting repeatedly for updating (by 
Eq. (0)) the largest w^s. For instance, the probability of 
obtaining a value r > KM (with K a small integer) is 
roughly the probability of extracting out of the r steps at 
least 2K times steps of size at least M/2. For K « r 
this is basically the probability of at least one step of size 
at least M/2 multiplied by r and risen at the power 2K: 



2 A' 



R{> r = KM,t) ~ [tS{> M/2)} 
~ [T(M/(2s min )) 



-al2K 



which can also be written as 



R{> r,r) 



N 



(2q)*l 



-r/M 



(27) 



i.e., R(> r, r) decays exponentially with r. Values larger 
than KM are still rigorously disallowed. 

To sum up the results for small and moderate intervals 
T < T max : except for the central dome region |r| < Tdome 
and the extremely large \r\ > M region, the probability 
distribution function i?(> r, r) is similar to the single step 
probability R{> r, 1) = S(> r) (Eq. @). This is intu- 
itively explained by the fact that the probability to arrive 
after r random steps at large (but less than M) r values 
is dominated by the probability of having a single step of 
order r. 

The 3 regions above: central dome region r < rdome, 
the Levy-like power law region, and the extremely large 
r > M exponential region are the main features of the 
distribution R(> r, r) for intervals r < r max . 



2.4 Cross-over for t > r max 

The time evolution of the R(> r, r) shape depends on the 
fact that while the upper cut-off region (beyond which 
the power law fails) r > M is fixed, the "dome" region 
expands with r according to Eq. (23). (Intuitively this is 



because, as one increases the number of time steps r, one 
can reach larger values of their sum r(r) even if each of 
the individual steps is of order w or less). As seen above, 



this leads for intervals r > T max (Eq. (g6|) ) to the disap- 
pearance of the intermediate Levy-like power law region 
in the R(> r, r) distribution. After this time, the central 
dome will keep expanding on the expense of the r > M 
cut-off region. 

As it expands, the dome will assume a shape closer 
and closer to a Gaussian. This will be consistent with the 
central limit theorem as the involved elementary steps will 
be ultimately distributed on a finite support s m j„ < s < 
M of quite limited extent compared with the r>M range 
of values probed by the dome for very large times. 

Indeed, for time intervals r 3> r max , the probability of 
many steps of size close to M is not negligible. Thus the 
dynamics consists in a random walk of individual steps of 
size distributed within the finite support between and 
M. The distribution becomes a Gaussian whose expansion 
is dominated by the largest steps M. Since there will be 
roughly one step of size M per r max = N a interval, the 
Gaussian width will expand as 



<j(t) = M{T/T rnax ) 1 l' 1 = MN~ 



a/2 T l/2_ 



(28) 



For returns much larger than this, i.e., r > ct(t), the ex- 
ponential regime Eq. ( p7|) will still survive. 

The above results can be verified by the numerical sim- 
ulations shown below. Moreover, for time interval r > 
T m ax but not too large, the power law distribution of re- 
turns with exponent well outside the stable Levy regime 
of a < 2 can be obtained in the simulations, with the 
exponential cut-off effect. 



3 Numerical simulations 

We have performed the computer simulations of autocat- 
alytic system described by Eqs. @ and ([?]), where the 
random factor A is set to uniformly distribute in the range 
0.9 < A < 1.1. In our simulations the number of subsys- 
tems N = 500, and the lower bound factor q = 0.3. Thus, 
according to relation (^), the exponent a of the power law 
distribution of w is about 1.4, which has been verified by 
previous simulations Jl2[ . 

Here we numerically study the distribution of the fluc- 
tuations or "returns" 



r(r) = lnW(t + r) -lnW(t) 



(29) 



for different time intervals r, to compare with the above 
analytic results. Note that in Eq. (g9|) we use the loga- 
rithmic difference for the definition of return, as in usual 
financial applications, which is approximately the relative 
change Eq. ( |l2] ) if the change is small. Note that due to 
the lower bound qui (Eq. ([?])), in this system W(t) has an 
increasing trend, which makes the return ( |29| ) more pos- 
sible to be positive, especially for large r. Thus, in our 
simulation results shown below, a maximum at positive 
finite r is obtained for large r distribution. In general, 
one could normalize W(t) (wi(t)) by a value exu(Kt/N) 
with k constant, for detrending, and then Eq. (129) would 
change by a constant: r(r) — > r(r) — kt/N, which does 
not influence the statistical properties. 



Z.F. Huang and S. Solomon: Power, Levy, Exponential and Gaussian Regimes 



5 




-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 

Returns r(t) 



Fig. 1. Semi-plot of the probability density of returns r(r) (defined by Eq. (|29|)) for different time intervals r = 1, 500, and 
f2000 (TV = 500). 




Returns r(t) 



Fig. 2. The probability density of the positive returns r(r) (defined by Eq. in log-log plot for time intervals r — 1, 50, 

200, and 500 (TV = 500, calculated after t = 10 5 TV updatings), with a straight line of slope —3 (a — 2) for comparison. 
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With N = 500 and a ~ 1.4, we have r max = N a ~ 
6000 from Eq. (prf), and the behaviors obtained for small 
and large r are different, as shown in Fig. |l| for r = 1, 500, 
and 12000. With the increase of time interval r, the sharp 
peak of the distribution curve is smeared out into a dome- 
like shapej similar to the empirical findings of financial 
markets 



3.1 T <C T max 

The numerical results of probability density R(r, r) for 
time intervals r = 1 (pluses), 50 (circles), 200 (diamonds) 
and 500 (stars) are presented in Fig. ||. The measurement 
was performed after t = 10 5 N updatings, and averaged 
over 100 runs for r = 1 and 800 runs for r > 50. The 
results for r 



1 are just what we derive in Sec. 2/2 



a straight line in log-log plot with sharp end, obeying Eq. 
( |l8| ) with a about 1.5, and a cut-off for large r. 

For intervals larg er than 1, the 3 regions obtained an- 
alytically in Sec. |7| are clearly shown in the log-log plots 
of Fig. g. The first one is the central "dome" -like region 
for small r, with larger extent for larger r, as predicted 
in Eq. (p3]). For small interval r = 50, the derivative at 
small r is not close to zero, that is, the sharp end persists, 
while it flattens into the dome for larger intervals of 200 
and 500. In the intermediate r range, the power law be- 
havior similar to Eq. (|l8|) is presented. For small values of 
t (< 200 in Fig. ||), the exponent a is within the stable 
Levy regime, that is, < a < 2, however, for larger r (say 
500) one can obtain the exponent a > 2, which is similar 
to the phenomenon shown below for r > T max . 

When the return r is large, the deviation from the 
straight line and the curvature in log-log plot can be ob- 
served (see Fig. ||), which is just the cut-off effect described 
above, that is, the exponential decay for far tail distribu- 
tion (similar to Eq. (|27|)). 



3-2 T > ^max 

The crucial results of this paper are for large time interval 
t, which corresponds to realistic time scale observed in 
nature (e.g., financial markets). These results are the con- 
sequences of the truncation in this autocatalytic system, 
as discussed analytically above, and can be compared with 
those of the real market data. 

We have calculated the return distribution for large 
r: 12000, 25000, and 50000 (averaged over 1000 runs), as 
shown in Fig. ||. Besides the dome-like shape for small r, 
what interests us is the power law region for intermediate 
and large returns. In this region the exponent a is about 
3, well beyond the stable Levy regime, but in agreement 
with the recent observations in real stock market HI). 

However, for extremely large returns, i.e., the far tail 
of the distribution, the exponential cut-off effect still re- 
mains, as shown in the bent of the log-log plots in Fig. |[ 
To see more clearly, we replot Fig. || on a semi-log scale 
to Fig. |[ where the tail of return distribution looks like 
a straight line, indicating an exponential-type behavior. 



This exponentially asymptotic decay was not detected in 
the empirical findings of H and 0, but agrees with the 
most recent observation in Hong Kong stock market as 
shown in Sec. || below. 

To see whether this phenomenon of a = 3 power law 
and exponential cut-off is related to finite size effect or 
is intrinsic, we perform the simulations for N = 50, and 
present the result of r = 2500 in Fig. || The range of power 
law behavior (a ~ 4.5 for N — 50) may be slightly shorter 
than that of the larger system TV = 500. Thus, the power 
law region before the exponential cut-off is expected to 
extend for large system size (N), which has been verified 
by simulations for N = 5000. 

For very large r, i.e., r 3> T max , the distribution of 
returns is to approach Gaussia n ba sed on the central limit 
theorem (also expected in Sec. 2.4), which has been found 
in reality ||. In Fig. || the result for r = 500000 is also 
shown with the Gaussian-like behavior for not too large r 
value, which exhibits as a parabola in semi-log plot (Fig. 
[I|) [|l6|. The exponential-type decay is still found in Fig. ^ 
for extremely large r as expected. 



4 Discussion and conclusion 

In finance, the tradings of the various investors Wi are per- 
formed independently. Therefore, the natural time mea- 
sure is not the number of operations but the number of 
operations divided by the number of components iV: 



T = t/N. 



(30) 



Thus, for the power law behavior of large t (500, 12000, 
25000, and 50000) shown in Figs. || and || with exponent 
well outside the Levy regime, the corresponding T values 
are 1, 24, 50, and 100 (for N = 500 in our simulations). 
Although in real market the time interval between trans- 
actions of stocks is irregular, one could estimate from the 
market transaction data [|17| that the unit scale T = 1 
here approximately corresponds to several minutes of real 
time. 

It is interesting to compare the results of this auto- 
catalytic system with actual measurements of the stock 
markets. The power law behavior with exponent a about 
3 was observed in recent empirical studies on S&P 500 H 
and German DAX [EJ, consistent with the results shown 
in Fig. U for intermediate and large returns, but not the 
exponential cut-off. However, very recently it has been 
found from the Hang Seng Index (HSI) of Hong Kong fit] 
that the index fluctuations for the first few minutes of 
daily opening behave very differently from those of the 
other times, due to much higher influences of exogenous 
factors at the opening. As shown in Fig. ||, if we skip the 
data in the first 20 minutes of each trading day, the dis- 
tribution for 1 minute time interval index move (defined 
as index(< + r) — indcx(i), with the statistical properties 
very similar to that of Eq. (E9T) for the high-frequency 
regime, e.g., interval r = 1 min) exhibits the phenomenon 
of exponential-type decay []l8| after the transient a = 3 
power law region, in agreement with our simulation results 
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Returns r(t) 



Fig. 3. The probability density of the positive returns r(r) in log-log plot for large time intervals r = 12000, 25000, 50000, 
and 500000 of N — 500, as well as r = 2500 of N = 50 (calculated after t = 10 N updatings), with a straight line of slope —4 
(a = 3) which fits the data for a few orders of magnitude. For r = 500000, only part of the distribution is shown. 




Returns r(t) 

Fig. 4. Replot of Fig. |^ on a semi-log scale, using larger bin size. 
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Fig. 5. Log- log plot of the probability density of 1 minute index moves for the Hong Kong HSI (from Jan. 1994 to Dec. 1997), 
with the skip of the first 20 minutes in daily opening J7J. A straight line with slope —4 is also shown. 



of Figs. U and ^ (r = 500). Moreover, real markets show 
the tendency of a crossover towards a Gaussian for long 
enough times |3|, which has also been found in our system 
(see Figs. || and || for very large interval r = 500000 and 
intermediate r). 

In order to account for the experimentally observed 
volatility correlations, one may follow Rcf. and feed 
back the absolute market returns (Eq. (^)) into the indi- 
vidual gain factor @. A possible form is using: 

A(i)=exp[(r)+r ? ((r 2 )-(r} 2 ) 1 / 2 ] 

where the angle brakets indicate averages over the last r 
steps and i] is a Gaussian random number of zero mean 
and unit variance. The numerical results for these more 
realistic simulations will be presented elsewhere. 

In summary, we have shown that the simple random 
multiplicative model of Ref. || , with a lower cut-off, gives 
many of the properties found in reality and in more com- 
plicated models, like e.g. the percolation model |l9|: power 
law with effective exponent a near 3, rounding of the sin- 
gularity at zero returns, crossover to Gaussian for long 
times. Wc still have to work on implementing volatility 
clustering, multifractality, lack of up-down and time re- 
versal symmetry, and correlation between traded volume 
and volatility. Only the original random walk model of 
Bachelier seems to us simpler than the present model but 
of the above properties it has only a Gaussian distribution. 
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